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4^ ! Abstract 

1 : P-tons produced in the early stage of non-central heavy-ion colhsions can develop a longitudinal 
I ^ i ' fluid shear because of unequal local number densities of participant target and projectile nucleons. 

I . Under such fluid shear, local parton pairs with non-vanishing impact parameter have finite local 

> : 

en , relative orbital angular momentum along the direction opposite to the reaction plane. Such finite 

g ; ... ... ... _ _ ... ... . ... ... 

• polarization along the same direction due to spin-orbital coupling. Local longitudinal fluid shear 

. is estimated within both Landau fireball and Bjorken scaling model of initial parton production. 

O ■ 

^ I Quark polarization through quark-quark scatterings with the exchange of a thermal gluon is cal- 



^ ' culated beyond small-angle scattering approximation in a quark-gluon plasma. The polarization 

is shown to have a non-monotonic dependence on the local relative orbital angular momentum 
dictated by the interplay between electric and magnetic interaction. It peaks at a value of relative 
orbital angular momentum which scales with the magnetic mass of the exchanged gluons. With 
the estimated small longitudinal fluid shear in semi-peripheral Au + Au collisions at the RHIC 
energy, the final quark polarization is found to be small \Pq\ < 0.04 in the weak coupling limit. 
Possible behavior of the quark polarization in the strong coupling limit and implications on the 
experimental detection of such global quark polarization at RHIC and LHC are also discussed. 



PACS numbers: 25.75.-q, 13.88.+e, 12.38.Mh, 25.75.Nq 



I. INTRODUCTION 



Collective phenomena and jet quenching as observed in high-energy heavy-ion collisions at 
the Relativistic Heavy-ion Collider (RHIC) at the Brookhaven National Laboratory (BNL) 
provide strong evidence for the formation of a strongly coupled quark gluon plasma 1, ^. 
Elliptic flow or azimuthal anisotropy of the hadron spectra in semi-peripheral heavy-ion 
collisions and its agreement with the ideal hydrodynamic calculations [3] indicate a near 
perfect fluid behavior of the produced dense matter. Such an empirical observation of 
small shear viscosity j4| is consistent with the large value of the jet transport parameter 
as extracted from jet quenching study of both single and dihadron spectra suppression ^. 
Study of the collective behavior is made possible by investigating hadron spectra in central 
rapidity region in non-central or semi-peripheral heavy-ion collisions. Extending the study 
to large rapidity region of non-central heavy-ion collisions should provide more information 
not only about the initial condition for the formation of the dense matter {oj] but also the 
dynamical properties of the strongly coupled quark-gluon plasma. 

Considering the longitudinal momentum distribution at various transverse positions in 
a non-central heavy-ion collision, one will find a longitudinal fluid shear distribution rep- 
resenting local relative orbital angular momentum. Recently, it has been pointed out that 
the presence of such local orbital angular momentum of the partonic system at the early 
stage of non-central heavy-ion collisions can lead to a global polarization of quarks and anti- 
quarks 0] in the direction orthogonal to the reaction plane. Understanding the spin-orbital 
interaction inside a strongly coupled system can open a new window to the properties of 
quark-gluon-plasma (QGP). Although no detailed calculations have been carried out, an 
estimate using a screened static potential model in the small angle approximation shows 
qualitatively that spin-orbital coupling in Quantum Chromodynamics (QCD) can lead to a 
finite global quark and anti-quark polarization. Such a global quark/anti-quark polarization 
should have many observable consequences such as global hyperoiipolarization [3, Q] and 
vector meson spin alignment 9]. Predictions have been made [3, 3 fo'^ these measurable 
quantities as functions of the global quark polarization Pg in various hadronization scenarios. 
Since the reaction plane in heavy-ion collisions can be determined in experiments by measur- 
ing the elliptic and direct flows, measurements of the global hyperon polarization or vector 
meson spin alignment become feasible. These measurements at RHIC are being carried out 



and some of the preliminary results have already been reported 

The estimate of the global quark polarization in Ref. ^ was obtained by evaluating the 
polarization cross section in the impact parameter space with small angle approximation in 
an effective potential model. The analytical result, 

Pg = -7Cfip/2E{E + m,) (1) 

has an intuitive expression, where p is the average cm. momentum of two partons with 
an average transverse separation due to the longitudinal fluid shear. However, for a 
massless quark in a small longitudinal fluid shear, the obtained quark polarization Pg can 
become larger than 1, indicating the breakdown of the small angle approximation. A more 
realistic estimate in non-central heavy-ion collisions at RHIC indicates a small value of the 
average longitudinal fluid shear. Therefore, it is imperative to have a more realistic estimate 
of the quark polarization Pg beyond the small angle approximation. This will be the focus 
of this paper. 

The rest of the paper is organized as follows. In Sec. II, we calculate the average longitu- 
dinal fluid shear in two different models of parton production. In a Landau fireball picture, a 
wounded nucleon model of bulk parton production in heavy-ion collisions is used with both 
simple hard-sphere and more realistic Wood-Saxon nuclear geometry. In the Bjorken scaling 
scenario, we use HIJING Monte Carlo model to estimate the transverse shear of the rapidity 
distribution of the produced parton in heavy-ion collisions at the RHIC energy which will be 
used to estimate the longitudinal fluid shear in the local comoving frame of the plasma. In 
Sec. HI, we use the Hard Thermal Loop (HTL) resummed gluon propagator in the comov- 
ing frame of the local longitudinal fluid cell to extend the calculation of quark polarization 
in Ref. [?! beyond small angle approximation and discuss the relative contributions from 
electric and magnetic part of quark-quark scattering. Finally in Sec. IV, we discuss the 
numerical results and their implications for experimental measurements at RHIC. 



II. ORBITAL ANGULAR MOMENTUM AND SHEAR FLOW 

Let us consider two colliding nuclei with the projectile of beam momentum per nucleon 
Pin moving in the direction of the illustrated in Fig. [H The impact parameter b, 

defined as the transverse distance of the center of the projectile nucleus from that of the 
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FIG. 1: Illustration of a non-central heavy- ion collision with impact parameter b. The global 
angular momentum of the produced dense matter is along —y, opposite to the reaction plane. 

target, is taken to be along the x-direction. The normal direction fib of the reaction plane, 
given by, 

nb = PinXb/\piny<b\, (2) 

is along y. For a non-central collision, the two colliding nuclei carry a finite global orbital 
angular momentum Ly along the direction orthogonal to the reaction plane {—y). How such 
a global orbital angular momentum is transferred to the final state particles depends on the 
equation of state (EOS) of the dense matter. At low energies, the final state is expected to be 
the normal nuclear matter with an EOS of rigid nuclei. A rotating compound nucleus can be 
formed when the colliding energy is comparable or smaller than the nuclear binding energy. 
The finite value of the total orbital angular momentum of the non-central collision at such 
low energies provides a useful tool for the study of the properties of superdeformed nuclei 
under such rotation 



17j. At high colliding energy at RHIC, the dense matter is expected 
to be partonic with an EOS of the quark-gluon plasma. Given such a soft EOS, the global 
orbital angular momentum would probably never lead to the global rotation of the dense 
matter. Instead, the total angular momentum will be distributed across the overlapped 
region of nuclear scattering and is manifested in the shear of the longitudinal flow leading 
to a finite value of local vorticity density. Under such longitudinal fluid shear, a pair of 
scattering partons will on average carry a finite value of relative orbital angular momentum 
in the opposite direction to the reaction plane as defined in Eq. ([2]). According to Ref. [3], 
quark (or antiquark) will acquire a global polarization after such scatterings through the 
spin-orbital coupling in QCD. 



The magnitude of the total orbital angular momentum Ly and the resulting longitudinal 
fluid shear can both be estimated within the wounded nucleon model of particle production 
in which the number of produced particles is assumed to be proportional to the number of 
participant nucleons. The transverse distributions (integrated over y) of participant nucleons 
in each nucleus can be written as, 

= j dydzpY{x, y, z, b), (3) 

in terms of the participant nucleon number density p^^{x, y, z, h) in nucleus A in the coordi- 
nate system deflned above. The superscript P or T denotes projectile or target respectively. 
The total orbital angular momentum Ly of the two colliding nuclei can be deflned as, 

Ly = -P^n j X dx{^^ d^^- 

where is the momentum per incident nucleon. 
We assume the Woods-Saxon nuclear distribution, 

fwsi^^ y^^M = C \ l + exp ^'^^ ^ ^^^^ ^ ^ — ^1 . (5) 



The participant nucleon number density is then 



PAwsi^^ = fwsi^^ y^z,h)\l- exp 



-aNN I dzf^s{x,y,z,b) 



(6) 



where a^^ ~ 42 mb is the total cross section of nucleon- nucleon scatterings at the RHIC 
energy, C is the normalization constant and a is the width parameter set to a = 0.54 fm. 

Shown in Fig. [2] as the solid line is the numerical value of Ly as a function of b for the 
Woods-Saxon nuclear distribution. As a comparison, we also plot as the dashed line the Ly 
distribution with a hard-sphere nuclear distribution which was used in Ref. [t]. With the 
hard-sphere nuclear distribution, the participant nucleon density is given by the overlapping 
area of two hard spheres. 



PAjisi^^ y^ = fA^si^^ y^ b)e{RA - ^/{x±b/2y + y^ + z^), (7) 

and 

OA 

f^%si^,y,z,b) = ^JiRA - ^/{xTb/2y + y' + z^), (8) 



where Ra = 1.1 2^4^/^ is the nuclear radius and A the atomic number. We note that there are 
significant differences between two nuclear geometry in the total orbital angular momentum 
Ly in the overlapped region of two colliding nuclei. In both cases, the total orbital angular 
momentum is huge and is of the order of 10^ at most impact parameters. 
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FIG. 2: (Color online) The total orbital angular momentum of the overlapping system in Au + Au 
collisions at the RHIC energy as a function of the impact parameter b. The solid and dashed curves 
are from the Woods-Saxon and hard-sphere distributions, respectively. 

Since RHIC data indicate the formation of a strongly coupled quark-gluon plasma jl|, 
we can assume that a partonic system is formed immediately following the initial collision 
and interactions among partons will lead to both transverse (in x-y plane) and longitudinal 
collective motion in the quark-gluon plasma (QGP). The total orbital angular momentum 
carried by the produced system will manifest in the longitudinal fiow shear or a finite value 
of the transverse (along x) gradient of the longitudinal fiow velocity. How the total angular 
momentum is distributed to the longitudinal fiow shear and the magnitude of the local 
relative orbital angular momentum depends on the parton production mechanism and their 
longitudinal momentum distributions. We consider two different scenarios in this paper: 
Landau fireball and Bjorken scaling model. 

By momentum conservation, the average initial collective longitudinal momentum at 
any given transverse position can be calculated as the total momentum difference between 
participating projectile and target nucleons. Since the total multiplicity in A + A collisions is 
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FIG. 3: The average longitudinal momentum distribution pz{x,b) in unit of po = ^/s/[2c{s)] as 
a function of x/{Ra — b/2) for different values of b/RA with the hard sphere (upper panel) and 
Woods-Saxon (lower panel) nuclear distributions. 



proportional to the number of participant nucleons [18| , we can make the same assumption 
for the produced partons with a proportionality constant c(s) at a given center of mass 
energy ^/s. 

In a Landau fireball model, we assume the produced partons thermalize quickly and have 
a common longitudinal flow velocity at a given transverse position of the overlapped region. 
The average collective longitudinal momentum per parton can be written as 

where po = ^/s/[2c{s)]. The distribution pz{x,b; ^/s) is an odd function in both x and b 
and therefore vanishes at x = or 6 = 0. In Fig. [21 Pz{x, b; ^/s) is plotted as a function of 



X at different impact parameters b. We see that Pz{x,b; ^/s) is a monotonically increasing 
function of x until tfie edge of the overlapped region \x ±b/2\ = beyond which it drops 
to zero (gradually for Woods-Saxon geometry). 

From pz{x,b; y/s) one can compute the transverse gradient of the average longitudinal 
collective momentum per parton dpz/dx which is an even function of x and vanishes at 
6 = 0. One can then estimate the longitudinal momentum difference Ap^ between two 
neighboring partons in QGP. On average, the relative orbital angular momentum for two 
colliding partons separated by Ax in the transverse direction is ly = —{AxYdpz/dx. With 
the hard sphere nuclear distribution, ly is proportional to dpo/dx = Pq/Ra = \/s/['^c{s)Ra]- 

In Au + Au collisions at ^/s = 200 GeV, the number of charged hadrons per participating 
nucleon is about 15 [18]. Assuming the number of partons per (meson dominated) hadron is 
about 2, we have c{s) ~ 45 (including neutral hadrons). Given Ra = 6.5 fm, dpo/dx ~ 0.34 
GeV/fm and we obtain Iq = —{Ax)'^dpo/dx ~ —1.7 for Ax = 1 fm. In Fig. HJ we show the 
average local orbital angular momentum ly for two neighboring partons separated by Ax = 1 
fm as a function of x for different impact parameter b for both Woods-Saxon and hard-sphere 
nuclear distributions. We see that ly is in general of the order of 1 and is comparable or 
larger than the spin of a quark. It is expected that c{s) should depend logarithmically on 
the colliding energy y/s, therefore ly should increases with growing y/s. 

In a 3-dimensional expanding system, there could be strong correlation between longi- 
tudinal flow velocity and spatial coordinate of the fluid cell. The most simplified picture 
is the Bjorken scaling scenario 1^ in which the longitudinal flow velocity is identical to 
the spatial velocity t] = log[(t + z)/(t — z)]. With such correlation, local interaction and 
thermalization requires that a parton only interacts with other partons in the same region 
of longitudinal momentum or rapidity y. The width of such region in rapidity is determined 
by the half-width of the thermal distribution f{Y,px) = exp[—pTCosh{Y — r])/T] [20], which 
is approximately Ay ^1.5 (with (pt) ~ 2T). The relevant measure of the local relative 
orbital angular momentum between two interacting partons is, therefore, the difference in 
parton rapidity distributions at transverse distance of (5x ~ l//x on the order of the average 
interaction range. 

One needs a dynamical model to estimate the local rapidity distributions of produced 



partons. For such a purpose, we use HIJING Monte Carlo model 
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221 ] to calculate the 



hadron rapidity distributions at different transverse coordinate (x) and assume that parton 




FIG. 4: The average orbital angular momentum ly = —{Ax)'^dpz/dx of two neighboring partons 
separated by Ax = 1 fm as a function of the scaled transverse coordinate x/{Ra — b/2) for different 
values of the impact parameter h/ Ra with the hard-sphere (upper panel) and Woods-Saxon (lower 
panel) nuclear distributions. 

distributions of the dense matter are proportional to the final hadron spectra. Shown in 
Fig. [5] is the average rapidity (Y) as a function of the transverse coordinate x for different 
values of the impact parameter h. The distributions have exactly the same features as 
given by the wounded nucleon model in Fig. [3l The variation of the rapidity distributions 
with respect to the transverse coordinate is illustrated in Fig. [6] by the normalized rapidity 
distributions 

at different transverse coordinates, x = 0, ±2 fm. At finite values of the transverse coordi- 
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FIG. 5: (Color online) The average rapidity (Y) of the final state 



transverse coordinate x from HIJING Mont Carlo simulations 
collisions at -v/i = 200 GeV. 
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FIG. 6: (Color online) The normalized rapidity distribution of particles P(Y,x) [Eq. (jlOp ] at 
different transverse position x from HIJING simulations of non-central Au + Au collisions at ^/s = 
200 GeV. 

nates x, the normalized rapidity distributions evidently peak at larger values of rapidity 
The shift in the shape of the rapidity distributions will provide the local longitudinal fluid 
shear or finite relative orbital angular momentum for two interacting partons in the local 
comoving frame at any given rapidity Y. To quantify such longitudinal fluid shear, one can 



calculate the average rapidity within an interval Ay at Y , 

+ (11) 

The average rapidity shear or the difference in average rapidity for two partons separated 
by a unit of transverse distance Ax = 1 fm is then, 

d{Y) ^ Al dHog P{Y,x) 
dx ^ 12 dYdx ■ ^ ^ 




FIG. 7: (Color online) The average rapidity shear d{Y)/dx within a window Ay = 1 as a function 
of the rapidity Y at different transverse position x from HIJING calculation of non-central Au + Au 
collisions at -y/i = 200 GeV. 

Shown in Fig. [7] is the average rapidity shear as a function of the rapidity y at different 
values of the transverse coordinate x for Ay = 1. As we can see, the average rapidity shear 
has a positive and finite value in the central rapidity region. The corresponding local relative 
longitudinal momentum shear is 

(13) 

ox ox 

With (pt) ~ 2T ~ 0.8 GeV, we have d{pz)/dx ~ 0.003 GeV/fm in the central rapidity region 
of a non-central Au + Au collision at the RHIC energy given by the HIJING simulations, 
which is much smaller than that from a Landau fireball model estimate. 



III. GLOBAL QUARK POLARIZATION 



As we have discussed earlier, under the longitudinal fluid shear, a pair of interacting 
parton will have a finite value of relative orbital angular momentum along the direction 
opposite to the reaction plane. In this section, we will calculate quark polarization via 
scatterings with fixed direction of the relative orbital angular momentum. We will assign a 
fixed direction of the impact parameter Xt between two interacting partons to refiect the 
direction of the relative orbital angular momentum. The magnitude of the relative orbital 
angular momentum will be charaterized by the relative longitudinal momentum p between 
two partons separated by a transverse distance Ax ~ 1/// on the order of the average 
interaction range. With the averaged longitudinal fluid shear dp^/dx in the center of mass 
frame of the two colliding nuclei in the Landau fireball model, we have p = Ax{dpz/dx). 
In the Bjorken scaling scenario with strong correlation between spatial and momentum 
rapidity, the average local longitudinal shear in the comoving frame will be given by p = 
AxpT cosh{Y)d{Y) / dx, where pr is the average transverse momentum. 

A. Quark scattering at fixed impact parameter 

We consider the scattering Ai) + g2(-P2, A2) QiiPa, A3) + g2(-P4, A4) of two quarks 

with different flavors, where Pi = {Ei^pi) and Aj in the brackets denote the four momenta 
and spins of the quarks respectively. The cross section in momentum space is given by, 

where Ai{Q) is the scattering amplitude in momentum space, Q = P3 — Pi = P2 — P4 is 
the four momentum transfer, Cqq — 2/9 is the color factor, and F — 4^/ (Pi • P2)^ — m\m\ is 
the flux factor. Since we are interested in the polarization of one of the quarks gi after the 
scattering, we therefore average over the spins of initial quarks and sum over the spin of the 
quark ^2 in the final state. 

We work in the center of mass frame of the two quark system. For simplification, we 
neglect thermal momentum in the transverse direction and assume the relative momentum of 
the two quarks separated by a transverse distance Ax of the order of the effective interaction 
range 1/// is simply given by the longitudinal fiuid shear. 



One can integrate over the final momentum p4 of tlie quark q2 and tlie longitudinal 
component psz of the quark qi, and obtain 



da 



As 



^qq 



4F4 



Ai,A2,A4 i=+,- 



-A<(g,)A^*(Q, 



32 I 



(15) 



where pf^ = -^^J — g^^, corresponding to two possible solutions of the energy-momentum 
conservation in the elastic scattering process, p = |pi | = | p2 1 ; and gr = Pst is the transverse 
momentum transfer. For simplicity, we will suppress the summation notation over z = +, — 
hereafter but keep in mind that the final cross section includes the two terms. 

Since we would like to calculate the polarization of one final-state quark with a fixed 
direction of the orbital angular momentum, or fixed direction of the impact parameter, we 
will cast the cross section in impact parameter space by making a two dimensional Fourier 
transformation in the transverse momentum transfer q^^ i.e.., 



d^a\^ 
d^XT 16F 



E 

Ai,A2,A4 



d^qr d'^kx 
(27r)2 (27r)2 



,i{kT-qT)-x: 



M{qT)M*{kT) 
Mot) A*{kT) ' 



(16) 



where A4{qT) and Ai{kT) are the scattering matrix elements in momentum space with four 
momentum transfer Q = (0,q) and K = (0, k) respectively, and 



A{qT) = ^JiE^ + E2)\pl\. 



(17) 



To calculate the quark-quark scattering amplitude in a thermal medium, we will use Hard 
Thermal Loop (HTL) resummed gluon propagator 23 



+ (a 



(18) 



where Q denotes the gluon four momentum and a is the gauge fixing parameter. The 
longitudinal and transverse projectors Pjt*^ are defined by 

-1 



(19) 
(20) 



with LJ = Q-U, = — U^u, q = —Q , g^u = Qi^w — U^U^. Here U is the fluid velocity 
of the local medium. The transverse and longitudinal self-energies are given by 23 1 



l--ln 
2 



1 + X 

1 — X 



TT 

+ ^-x 



y + ^(l-x2)ln 



X 



X 



X 



TT 



i—x(l — X 
4 ^ 



(21) 
(22) 



where x = uo/q and = g'^{Nc + A^//2)T^/3 is the Debye screening mass. 

With the above HTL gluon propagator, the quark-quark scattering amphtudes can be 
expressed as 

Miqr) = uxsiPi + Qh^ux,iPi)A^''{Q)uxAP2 - Q)luUx,{P2), (23) 

M*{kT) = ux,{Pihc.ux,{Pi + K)A"^*{K)ux,{P2h/3UxAP2 - K). (24) 
The product M.{qT)M.*{kT) can be converted to the following trace form, 

M{qT)M\kT) = A'^^(Q)A"'^*(ir)TrK(Pi + K)ux,{Pi + Q)l,{fi+ "^i)7a] 

Ai,A2 

xTrK(i'2 - K)uxAP2 - Q)7.(^2+ ^2)7/3]- (25) 



In calculations of transport coefficients such as jet energy loss parameter 24| and thermal- 



ization time [23] which generally involve cross sections weighted with transverse momentum 
transfer, the imaginary part of the HTL propagator in the magnetic sector is enough to 
regularize the infrared behavior of the transport cross sections. However, in our following 
calculation of quark polarization, total parton scattering cross section is involved. The con- 
tribution from the magnetic part of the interaction has therefore infrared divergence which 



can only be regularized throug 



1 the introduction of non-perturbative magnetic screening 



mass /i^ ^ Q.2hh^Nj2g'^T |26|. 

Since we have neglected the thermal momentum perpendicular to the longitudinal flow, 
the energy transfer = in the center of mass frame of the two colliding partons. This 
corresponds to setting x = in the HTL resummed gluon propagator in Eq. (fT8|) . In this 
case, the center of mass frame of scattering quarks coincides with the local comoving frame 
of QGP and the fluid velocity is U'^ = (1,0,0,0). The corresponding HTL effective gluon 
propagator in Feynman gauge that contributes to the scattering amplitudes is reduced to, 

\ 

The differential cross section can in general be decomposed into a spin-independent and 
a spin-dependent part, 

^^^A3 _ _do_ ^ ^ dAo_ 
(Pxt (Pxt (Pxt 

with 



A^-(g) = \2^,. (26) 



d^XT 2 yd'^XT d^XT J 



(PAct 1 / dcr^ da_ \ , . 

The spin-dependent part will mostly determine the polarization of the final state quark 
qi via the scattering. The calculation is involved. A simple estimate was given in Ref. 
using a screened static potential model and small angle approximation. In this case, the 
cross sections can be written in a general form as, 

d'^a , , , , 

F{XT,E), (30) 



d'^XT 

^ = n-ixTy<p)AF{xT,E), (31) 

where n is the polarization vector for qi in its rest frame. F{xt, E) and AF{xti E) are 
functions of both xt = \xt\ and the cm. energy E of the two quarks. We can show that 
the quark-quark scattering with HTL propagators has the same form as that in the static 
potential model But the detailed expressions of F{xt, s) and AF(xt, s) are much more 
complicated. 

In fact, one can show that these two parts of the cross sections should have the same form 
as given in Eqs. (|30l) and ( ISTl) due to parity conservation in the scattering process. We note 
that in an unpolarized reaction, the cross section should be independent of any transverse 
direction. Hence da/d'^XT depends only on the magnitude of xt but not on the direction. 
For the spin-dependent part, the only scalar that we can construct from the available vectors 
is n ■ (p X Xt)- 

We note that, x p is nothing but the relative orbital angular momentum of the two- 
quark system, I = xt 'x P- Therefore, the polarized cross section takes its maximum when n 
is parallel or antiparallel to the relative orbital angular momentum, depending on whether 
AF is positive or negative. This corresponds to quark polarization in the direction / or — /. 

As discussed in the last section, the average relative orbital angular momentum I of two 
scattering quarks is in the opposite direction of the reaction plane in non-central A + A 
collisions. Since a given direction of I corresponds to a given direction of xt-, there should 
be a preferred direction of xt at a given direction of the nucleus-nucleus impact parameter 
h. The distribution of xt at given h depends on the collective longitudinal momentum 
distribution shown in the last section. For simplicity, we consider a uniform distribution of 
Xt in all possible directions in the upper half xy-plane with x > 0. In this case, we need 
to integrate dAa /d'^XT and da/d'^XT over the half plane above ?/-axis to obtain the average 



cross section at a given b, i.e. 



= X '"X. * '''' 



The polarization of the quark is then obtained as, 



P, = ^. (34) 



B. Small angle approximation 

We only consider light quarks and neglect their masses. Carrying out the traces in 
Eq. (l25l) . we can obtain the expression of the cross section with HTL gluon propagators. 
The results are much more complicated than those as obtained in Ref.p] using a static 
potential model. However, if we use small angle or small transverse momentum transfer 
approximation, the results are still very simple. In this case, with ~ and = l^rl 
we obtain the spin-independent (unpolarized) cross section , 

<P<y g^Cqq f (Pqr Pkr 1 , 1 1 , 1 



+ „2 , ,.2 )( l2 , ..2 + l2 , ..2 )' (3^) 



and the spin-dependent differential (polarized) cross section. 



i2 



d'^XT 16 J (27r)2 (27r)2 p 

x(^ ^ + - t)(t2 2" + 1:2 r)- (36) 

We note that the polarized differential cross can be related to the unpolarized one by, 

xn).V.^. (37) 
Completing the integration over the transverse momentum transfer. 



2 ' 



(27r)2g2+^2^ J 271 

where Jo and Kq are the Bessel and modified Bessel functions respectively, we obtain, 

d'^a c, 



(38) 



qrdqT , 2 = KoiUmXr), (39) 



d^XT 2 



^a,^[Ko(/in^^T) + Ko{fiDXT)]% (40) 



[Ko{fl^XT) + Ko{flDXT)][fJ^mKi{fimXT) + /^D^l (yUDa^r)] , (41) 



d'^XT 2 p"^ 

where xt = xt/xt is the unit vector of xt- We compare the above results with that in the 
screened static potential model (SPM) where one also made the small angle approximation, 



da 



d^XT 



g Ct f d qx d kx 



,i{kT-qT)-XT _ 



(42) 



SPM 8 J (27r)2(27r)2 ^2(^2 + ^2^)(^2 + ^2^) • I ) 



SPM 4 J (27r)2(27r)2 g2 + ^2^ + /^s^ ' 

-d^XT- 

We see that the only difference between the two results is the additional contributions from 
magnetic gluons, whose contributions are absent in the static potential model. Using Eqs. 
(l38ll and (l39ll . we recover the results in Ref. [3], 

da 



d^XT 



SPM 



a^crK^lfiDXT), (44) 



dAa 



(pxn)- Xt 



O^sCt 5 fiDKo{flDXT)Ki{flDXT))- (45) 



-d'^XT^SPM '"^ p^ 

C. Beyond small angle approximation 

Now we present the complete results for the cross-section in impact parameter space using 
HTL gluon propagators without small angle approximation. The unpolarized and polarized 
cross section can be expressed in general as, 



dcT _ g Cqq f d qr d kx i,kT-qT) ST f^^T.h 



d^XT 64p2 J (27r)2 (27r)2 A{qT)A{kT) 
dAa g\g f d'^qr d'^kr Mk^^g^y^^ Af{qT,kT) 



(46) 



^i(kT~qT)-XT _i_}2±±_±_!_ ^47"^ 



d^XT 64p3 J (27r)2 (27r)2 A(gT)A(fcT) 



where the kinematic factor becomes A(g7^) = A/2p[p^; / and Af are given by, 

L \ _ Amm{qT,kT) A^eiqT.kx) 



Ame{qT,kT) , Ame{kT,qT) 



(q' + filW + fil) {q' + filW + fil 



.r,^ AAmm{qT,kT) , AAee{qT,kT 

Af{qT,kT) = 7^ — o I — ^ + 



(48) 



{q' + f^lW + f^D {q' + f^lW + t^D 

AAmeiqT,kT) AA^e(^T, ^t) 



{q2 + ^l){k^ + ^l) (g2 + ^2^)(p + ^2^ 



(49) 



A^UqT,kT) = {q-kf + 9,p\q-k) + Sp\q, + k,) + lQp^; (50) 

Aee{qT. kr) = Ammiqr, kr) + MQz + k,)[iQ- k) + p{q, + k,) + 2p^]- (51) 
Ame{qT, kr) = Amm{qT, kr) + [-2q^k^{q- k) + Apk^{q- k) - 2pqlk^ 

-2pq.,k,^ + ApV + ^P^k,)]] (52) 
/\Amni{qT,kT) = -{[q ■ k + Ap"^ - 2p{q^ + k^)]{k^qT - qzkr) 

+2p{q-k + Ap'^){qT-kT)}-{pxn); (53) 
Av4ee(gT, kr) = AAmm{qT, ^t) " 4p(g^ + kz)[{kzqT - qzkr) - p{kT - qr)] ■ (p x n); (54) 
AAmeiqr, kr) = AAmm{qT, kr) + 2pkz[2p{kT - qx) + {qz - kz)qT 

-{kzqr - qzkr)] ■ (p x n) + 2qzkz{kzqT - qzkr) ■ (P x n), (55) 

where p = \p\. It is useful to note that 

Ammiqr, kr) = A^^Q^t, qr), (56) 

Aeeiqr, kr) = A^eikr, qx)- (57) 
Hence, /(gr, kx) is symmetric in its two variables 

f{qxM) = fikx.qx). (58) 

Similarly from 

AAmm{qx, kr) = -AAmmikr, qx), (59) 
AAeeiqx, kx) = -AAeeikr, qx), (60) 
we know that Af{qT, kx) is anti-symmetric, 

Af{qT,kx) = -Af{kT,qx). (61) 

As mentioned above, to get the average polarization for a fixed direction of the reaction 
plane in heavy-ion collisions, we need to average over the distribution of xr- For this 

n 

purpose, we take the approach as in Ref. [7|], and integrate da /(Pxt and dAa /(Pxt over the 
half plane above the y-axis as shown in Eqs. (132!) and (133!) . It is convenient to carry out 
first the integration over x and y then that over qx and kx- To do this, we use the identity, 



where V denotes the principal value. 

It is useful to note that 5'^ikT — QT) = S'^iQT — kT) and V-rr^ — = —V — Therefore, only 

f^x Qx Qx f^x 

the first term on the r.h.s. of Eq. (1621) contributes to the total unpolarized cross section, 

^ g^Cqq 1 r (Pgr figT^qr) ^ . 

^ 6V2 4<p(2vr)2 A2(g^) " ^^^^ 

The polarized cross section Aa receives contribution only from the second term. 



Aa = ^ !^^' ^ /-V^^ dk, 1 Af{g,,gy;k,,gy) 

64p3 27r J_^^-2Z-^ 2ti J.^/^iZ^ 2ti k^ - g^ A(gr)A(A;r) 

Changing the integration variable gx = psin ^ and ^ = sin2(^/2) in the expression of the 
total cross section a, we obtain, 

a^^ r,J^±iU + ^l^ + . j, (65) 

where = {^lD/Tf = 47r«,(Ar, + Ar^/2)/3, Pm = {f^jTY = 0.255\A7ryalNj2, and 
T = T/^/§ with the center of mass energy of the gg-system. The integration can be 
carried out analytically, 

a ^^{4 + ^ + ^ ^-2(/3.T^ + ^±M!^)ln(l + ^) 

_2(/?^f 2 + PmT'-2PnT[-l 

Similarly, we make the variable substitutions g^ = — t^, g^ = pt\/l — /c^; = 



a/1 — in the integration for Aa and obtain, 

c.."? r r ^^v^ 



Aa = — ^ / dt dri 



X 



8vrs j_i Jo Jo vr^vT^vr^ 

;i - t')(4 + + tr^) - 2{t^^r] + 1) + 2t(l + ^r/)(5 + t'e^)/(e + ^7) 
il-t^ + 2f3mf^)il-tr] + 2P^f^) 

^ (1 - t^)(te + tr]) + 2(t2e^ + 1) + 2t{i + er7)(i + t^e^)/(e + v) 

il-t^ + 2Pnf^)il-tr] + 2PDf^) 
^ 2(1 - t^)(2 + + tr^) + 8t(l + er7)(l + tr7)/(e + v) \ .g^s 
(l-te + 2/?^r2)(l-t7^ + 2/3z,T2) /• 

Note that in the calculation of both the polarized and unpolarized cross sections we have 
limited the range of integration over the transverse momentum due to energy conservation. 



Such restriction is not imposed in the small angle approximation in Ref. 



3]. We see that, 



at given (3^ and (3d, both a and Aa are functions of the variable T = T/\/§. Since (3^ and 
Pd depend on as, the polarization Pg = Aa/cr also depends on the value of as- 

We can now carry out the integration numerically to get the quark polarization Pg. Before 
we show the numerical results, it is useful to look at two limits. 

(1) High energy limit. At very high energies, ^ T or T ^ 1, we have, 

a = + L_ + _i_to/Mi±ift^^. ,68) 

;i + 4&n& /3d -/3„ 13,41 + iHoT^) ' 




Act = - " I dt — t,Rn-\/ — : — - tRTi-^ ' 



(69) 

This is the case where the small angle approximation can be made. The above can also be 
obtained from Eqs. fl35tl36p given in the last section by carrying out the integration over xt 
in the half plane of x > 0. 

(2) Low energy limit. In the limit ^ T, we have + ~ and + /i^ ~ /i^, 
the cross sections become 

CqqOl ,\fl^A / 8 2 



8S "'^Hsa'^WO;'' 



Act = 



V192 2 '47 V192 '4' M 



+ f — ^^(-) + ^^(- 

Given the corresponding values of the F-function, one can obtain numerically, 

CqgaLv^x4 f 2.03 , 3.15 , 6.30 



(71) 



167rs T \ (3^ (3fj PmPo J 
We see that in the low energy limit the magnetic part contributes with different sign from 
the electric one. The polarization Pg = Aa/a is given by 

3 -2.03/?!, + 3.15/?^ + 6.30/?^/?o 

which tends to be a constant in this low energy limit. In the weak coupling limit <^ 1, 
/3d ^ I3m the above constant Pg 0.04 is a small positive number. 

It is also interesting to look at the contributions from the electric part only. The corre- 
sponding cross sections, denoted with subscription E, are 

da \ g\g f (fqr (Pkr A,,{qT,kT)e'^^^-'^^y^^ 



(PxtJe Q'^p' J (27r)2(27r)2A(g^)A(^^)(g2 + ^2^)(^2 + ^2^) 



(74) 



XT 



^ \Ap^ J (27r)2(27r)2A(g-'^)A(fc^)(g2 + ^2^)(p + ^2^)- 
Carrying out the integration over (Pxt in the half plane with x > 0, we obtain, 



As 



2 + ^ - 2(1 + PdT^) ln(l + ^) 



(75) 



(76) 



Aaf 



CqqOi 



X 



dt df dri — ; , , 

(1 - f'M + tv) + 2(t'er/ + 1) + 2t(l + ^r/)(l + t^^v)/{^ + v) 



(1 - + 2/5dT2)(1 -tr] + 2PdT^) 
In the high energy limit, where small angle approximation is applicable, we have. 



(77) 



Act, 



TlCqqOig 



4s/?dT2(1 + 4/5z,T2)' 
dt 

L t2 + 4/?^f 2 



4c /"^ 



tan 



t2 + 4/3,5T2 



In the low energy limit, we have. 



(Te 



AaE 



CqqOtg 



levrs T ' PI \192 M 



The polarization in this case. 



47r2 V192 



-0.24 



(78) 
(79) 



(80) 
(81) 



(82) 



is a negative constant which can also be obtained from Eq. (1731) by taking the limit Pm ^ Pd- 



D. Numerical results 



We now carry out the integration in Eq. fl67|) numerically and obtain the results for the 
quark polarization at intermediate energies between the high-energy and low-energy limit. 
The results are shown in Fig. Oas functions of V§/T. The quark polarization {—Pq) along 
the reaction plane approaches a small negative value as we have shown in the last subsection 
in the low-energy limit. The value of the low energy limit varies with as as given by Eq. 
( 173|) . Such a dependence on is a consequence of the magnetic and electric screening 
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FIG. 8: (Color online) Quark polarization -Pq as a function of \f%lT for different a^'s. 

masses in the polarized and unpolarized cross sections which have different dependence 
on Ug- However, from Eq. (173 p . the low-energy limit of the quark polarization becomes 
independent of in the weak coupling limit — > when /3m -C f^B- 

As one increases the relative cm. energy, the quark polarization changes drastically with 
VI/T. It increases to some maximum values and then decreases with the growing energy, 
approaching the result of small angle approximation in the high-energy limit. This structure 
is caused by the interpolation between the high-energy and low-energy behavior dominated 
by the magnetic part of the interaction in the weak coupling limit < 1. Therefore, the 
position of the maxima in ^ should approximately scale with the magnetic mass /^m- This 
is indeed the shown in Fig. [9l 

To further understand the interpolation between the high and low-energy limits in the 
numerical results, we also compare them in Fig. [10] with the results with the electric gluon 
exchange only. Without the contribution from the magnetic gluon interaction, the quark 
polarization takes a relatively large value Pq ~ —0.24 at low energies and then decreases 
with at high energies. The magnetic interaction in the low-energy limit apparently has a 
different sign in the contribution to the polarized cross section relative to that of the electric 
one. The net polarization is therefore reduced at finite to smaller negative values when 

<^ 1. The electric contribution to the net quark polarization also corresponds to the limit 
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FIG. 9: (Color online) Quark polarization -Pq as a function of \/I//im for different values of 



/im ^ yUz) or 1 in the full result. Even though our perturbative approach is no longer 

valid in such a limit, it indicates that the net quark polarization remains a finite negative 
value in the strong coupling limit as shown in Fig. [9l 

In Fig. [11] we also compare the full numerical results (solid lines) with those of the small 
angle approximation in the high-energy limit (dashed lines) as given by Eqs. fl68l) and fl69l) . 
These two groups of results indeed agree with each other at high energies. However, they 
both are different from the results of the static potential model in the small angle limit 
(dotted lines) 9l] which does not have the energy conservation restriction in the integration 
over the transverse momentum transfer. 



IV. CONCLUSIONS AND DISCUSSIONS 

In this paper, we have extended an earlier study 9| of the global quark polarization 
caused by the longitudinal fluid shear in non-central heavy-ion collisions. We have calculated 
the average local relative orbital angular momentum or longitudinal fluid shear with two 
extreme models: Landau fireball and Bjorken scaling model. In the Landau fireball model, 
we assumed a wounded nucleon model for local particle production with both the hard- 
sphere and Woods-Saxon nuclear distributions. Each parton is then assumed to carry an 
average longitudinal flow velocity calculated from the net longitudinal momentum at a given 
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FIG. 10: (Color online) Quark polarization -Pq as a function of ^f§/T with the full HTL gluon 
propagator (solid) as compared to the results with the electric part of the interaction only (dashed) . 
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FIG. 11: (Color online) Comparison of the results obtained using HTL gluon propagator (solid 
line) with those under small angle approximation (dashed line) and those using the screened static 
potential model under small angle approximation (dotted line). 



transverse position. In the Bjorken scaling model we considered correlation between spatial 
and momentum rapidity in a 3-dimensional expanding system for which we calculated the 
average rapidity or longitudinal momentum shear (derivative of the average rapidity or the 
longitudinal momentum) with respect to the transverse position x. The shear determines the 
local relative orbital angular momentum in the comoving frame at a given rapidity. These 
two model calculations provide estimates of the local fluid shear in two extreme limits. 

We have also extended the calculation of the global quark polarization Pg within per- 
turbative QCD at finite temperature beyond the small angle approximation of the previous 
study which might not be valid for small values of the local longitudinal fluid shear or 
the average cm. energy of a colliding quark pair. We found that the magnetic part 
of the interaction in one-gluon exchange is particularly important at low energies which 
cancels the contribution from the electric interaction and leads to smaller negative values 
of the net quark polarization in the weak coupling limit (a^ < 1). The final global quark 
polarization therefore is small in both the low and high-energy limits. It can, however, reach 
a peak value of about Pg ^ —0.24 at an energy determined by the nonperturbative magnetic 
mass V§ ~ 4/im, ~ g^T^y Nc/2. For \/§ < jj^m, the average quark polarization becomes 
significantly smaller. 

In semi-peripheral Au + Au collisions {b = Ra) at the RHIC energy, one can assume an 
average temperature T k, 400 MeV ^22]. With ~ 0.3, the global quark polarization reaches 
its peak value at cm. energy about 1.8 GeV. Since the magnetic interaction dominates the 
quark-quark interaction in our calculation, we can assume that the average interaction range 
in the transverse direction is given by the magnetic mass. Ax ~ 1/yUm- According to our 
estimates of the longitudinal fluid shear, the average cm. energy of the quark pair under 
such fluid shear is ~ 0.8 GeVV/im ~ 0.17 GeV (from Fig. ^ in the Landau fireball 
model. In the Bjorken scaling model (from Fig. [7]), the cm. energy provided by the local 
fluid shear is \fl ~ 0.004(pr) cosh(t/) imT^ / ^ 0.001 GeV in the central rapidity region 
(we assume (pt) ~ 2T). In both cases, the longitudinal fluid shear is so weak that the global 
quark polarization due to perturbative quark-quark scatterings is quite small according to 
our numerical calculations that go beyond the small angle approximation. 

In heavy- ion collisions at the Large Hadron Collider (LHC) energy ^/s = 5.5 TeV, the 
average multiplicity density per participant nucleon pair was estimated to be about a factor 



of 3 larger than that at the RHIC energy 28|]. The corresponding longitudinal fluid shear 



and the average cm. energy of a quark pair will be about a factor 6 larger than that at the 
RHIC energy in the Landau fireball model, assuming the temperature is about 1.44 higher. 
One can also expect the average local longitudinal fluid shear in the Bjorken scenario at 
LHC is similarly amplified compared to the RHIC energy in particular at large rapidity. 
Therefore, the resulting net quark polarization should also be larger at LHC. 

We want to emphasize that the above numerical estimate is based on a perturbative 
calculation via quark-quark scatterings in the weak coupling limit. It is still possible that 
quarks could acquire large global polarization through interaction in the strong coupling 
limit, as hinted by our results with large values of the strong coupling constant even though 
such a perturbative approach becomes invalid. The finite value of the quark polarization 
could be detected via measurements of the global hypcron polarization or the vector meson 
spin alignment with respect to the reaction plane. According to our estimate of the longitu- 
dinal fluid shear, the effect is more signiflcant at large rapidity under the Bjorken scenario 
of the initial parton production. 

In the limit of vanishing local orbital angular momentum provided by the longitudinal 
fluid shear, the approach we used in this paper in the impact-parameter representation 
might not be valid anymore. However, the final spin-polarization due to the spin-orbital 
interaction should approach to zero in this limit, which is approximately the result of our 
full calculation. 
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